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Abstract. The paper studies the question of whether the classical mirror and synchronous cou- 
plings of two Brownian motions minimise and maximise, respectively, the coupling time of the 
corresponding geometric Brownian motions. We establish a characterisation of the optimality of 
the two couplings over any finite time horizon and show that, unlike in the case of Brownian 
motion, the optimality fails in general even if the geometric Brownian motions are martingales. 
On the other hand, we prove that in the cases of the ergodic average and the infinite time horizon 
.^f criteria, the mirror coupling and the synchronous coupling are always optimal for general (possibly 

non-martingale) geometric Brownian motions. We show that the two couplings are efficient if and 
only if they are optimal over a finite time horizon and give a conjectural answer for the efficient 
couplings when they are suboptimal. 
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r^ ' 1. Introduction 

-i— > | 

Let the process B = (Bt)t>o be a fixed standard Brownian motion and consider a standard 
Brownian motion V = (Vt)t>o on the same probability space. For any starting points i,i/£l, 
define the coupling time t(V) to be the first time the processes x + B and y + V meet. It is obvious 
that the synchronous coupling V = B maximises the coupling time as it makes it infinite almost 
surely (assuming x 7^ y). Note further that the coupling time t(V) for any Brownian motion V 
cannot be smaller than the first time one of the processes x + B and y + V reaches level (x + y)/2. 

In the case of the mirror coupling V = —B, this random time actually equals t(V). Hence the 

O ■ 

mirror coupling minimises the coupling time of x + B and y + V almost surely. In particular, for 

any fixed T > 0, the extremal Brownian motion in the optimisation problem, 

(1) minimise (resp. maximise) P(r(V) > T) over all Brownian motions V, 

is given by the mirror (resp. synchronous) coupling, uniformly over all finite time horizons. 

It is natural to investigate the following closely related problem for geometric Brownian mo- 
tion: minimise the coupling time of the processes dXt = cr\Xt dBt and dYj(V) = o~2Yt(V) dVt over 
all Brownian motions V on a given filtered probability space. The aim here is to maximise the 
probability of the event that X and Y(V) couple before a given fixed time T. Since the processes 
X and Y(V) are, at any time t, given by the explicit deterministic functions of Bt and Vt respec- 
tively, the discussion above might suggest that mirror coupling of B and V should be optimal. 
Furthermore, since X and Y(V) are martingales, the Dambis-Dubins-Schwarz representation of 
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the difference X — Y(V) intuitively suggests that the two processes will meet as early as possible 
if the coupling of the Brownian motions B and V is chosen so that the instantaneous volatility of 
X — Y(V) is as large as possible. Equivalently put, the minimal coupling time should be achieved 
by the Brownian motion V which maximises (at every moment in time) the instantaneous qua- 
dratic variation d[X - Y(V)] S = ((^Xs) 2 + (a 2 Y s (V)f) ds - 2a 1 X s a 2 Y s {V) d[B, V] s . Since the 
(random) Lebesgue density of the covariation measure d[5, V] s on [0,oo) is always between — 1 
and 1, it follows that the mirror coupling V = —B should be optimal. However, as we shall see, 
both of these intuitive arguments turn out to be false in general. 

This paper investigates the problems of minimising and maximising the coupling time of two 
general (i.e. possibly non-martingale) geometric Brownian motions (GBMs) using a finite time, 
infinite time and ergodic average criteria. In the finite time horizon case we study the analogue 
of Problem ((T} for GBMs and give a necessary and sufficient condition on the value function for 
the mirror (resp. synchronous) coupling to be optimal. This leads to an if-and-only-if condition 
on the parameters of the GBMs, which characterises the suboptimality (and hence optimality) of 
the mirror (resp. synchronous) coupling for any finite time horizon. In contrast to the intuitive 
arguments given above, this condition implies that mirror (resp. synchronous) coupling can be 
suboptimal in Problem (pQ) for GBMs even if the geometric Brownian motions are martingales. 
This raises a natural question: is the exponential tail of the mirror (resp. synchronous) coupling 
optimal or, put differently, is the coupling efficient in the sense of [3]? We show that the mirror 
(resp. synchronous) coupling is efficient if and only if it is optimal, and hence may be inefficient. 
In the case where the coupling is suboptimal, the proof of the aforementioned equivalence suggests 
the conjecture that the synchronous (resp. mirror) coupling is efficient in the minimisation (resp. 
maximisation) problem. 

The stationary and infinite time horizon (for some "discount" rate q > 0) problems are given 
as the analogues of Problem (JT]) with P(r(V) > T) replaced by 

-i pT /*oo 

limsup- / F(T(V)>t)dt and / e" 9t P(T(y) > t) dt, 

respectively It is clear that in the case of Brownian motion, the mirror (resp. synchronous) 
coupling is optimal according to both of these criteria. In this paper we prove that, unlike in 
the finite time horizon case, the same holds for all (possibly non-martingale) geometric Brownian 
motions. In particular this implies that the mirror coupling, which may be inefficient (i.e. has 
a thicker exponential tail than the optimal coupling), nevertheless minimises both the Laplace 
transform of the tail probability for any "discount" rate q and its ergodic average. Our proofs are 
based on Bellman's principle. 

The mirror coupling and the synchronous coupling of Brownian motions and related processes 
have attracted much attention in the literature. For example the classical book [7] and paper [8] 
introduce the mirror couplings of Brownian motions and diffusion processes (see also book |10| 
for the general theory of coupling) . In [5] it is established that the mirror coupling is not the only 
maximal coupling, although it is the unique maximal coupling in the family of Markovian (also 
known as immersed) couplings. In [2] it is proved that the tracking error of two driftless diffusions is 
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minimised by the synchronous coupling of the driving Brownian motions. In [6j generalised mirror 
coupling and generalised synchronous coupling of Brownian motions are introduced; the former 
minimises the coupling time and maximises the tracking error of two regime-switching martingales, 
whereas the latter does the opposite. Articles p], [3], and [9] discuss various applications of the 
mirror coupling of reflected Brownian motions and other processes. In particular in |3J, the notion 
of efficiency of a Markovian coupling, also used in the present paper, is studied in the context of 
the spectral gap of the generator of a Markov process. 

The remainder of the paper is organised as follows. Section [2] describes the setting and basic 
notation, which is used throughout. Section[3]establishes the optimality of the mirror and synchro- 
nous couplings in the infinite time horizon (Section 13.11 Theorem [1]) and stationary (Section 13.21 
Proposition (5|) problems. In Section [J] we characterise the optimality of the mirror and synchro- 
nous couplings over a finite time horizon (Section 14. 1\ Theorem [SJ and analyse the efficiency of 
the two couplings (Section 14.21 Proposition [9|) . Appendix [A] contains a well-known lemma from 
stochastic analysis, which enables us to apply Bellman's principle. 

2. Setting and notation 

Let (Q, F, (J-'t)t>o,^') be a filtered probability space which is rich enough to support a standard 
(J-t)-Brownian motion B = (Bt)t>o- Let 

(2) V := {V = (Vj)t>o; V is an (J r i )-Brownian motion with Vq = 0} 

be the set of all (standard) (.F^-Brownian motions on this probability space. 

Let X = (X t )t>o and Y(V) = (Y t (y))t>o be geometric Brownian motions, satisfying stochastic 
differential equations 

(3) X t = x + [ X s (<ti dB s + a x da) and Y t (V) = y + [ Y S (V) (<r 2 dV s + a 2 ds) . 

Jo Jo 

The Brownian motion B is fixed throughout and V is any element of the set V, defined in ([2]). 
We assume throughout the paper that 

(4) x,y>0, 01,02 Gl and <Ti,o"2€M, such that o\02 > 0, 
and define the following constants 

(5) fM := a2 — a\ + o"i/2 — o\j1 and a± := 02 ± o\. 

Note that (gj) implies \a + \ > |cr_| > 0. 

Define the coupling time of the two processes in ([3]) as 

t(V) := inf{t > 0; X t = Y t (V)} (inf := oo). 

The random variable t(V) is zero when the two processes start at the same point and positive 
P-a.s. otherwise. Under mild assumptions (e.g. if the filtration {J-t)t>o is right-continuous or if all 
the paths of Brownian motions V and B are continuous), r(V) is P-a.s. equal to an (J r t )-stopping 
time. 
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3. Stationary and infinite time horizon problems 

3.1. Infinite time horizon problems. For any q > 0, we consider the following two problems: 
find V inf G V and F sup G V (if they exist) such that 

(qlnf) inf / e - 9t P(r(V) > t) dt = / e~ qt P(r(y inf ) > t) dt 

v&v Jo Jo 

and 

(qSup) sup/ e- qt F(T(V)>t)dt= e- qt ¥(T(V sup ) > t) dt. 

vgvJo Jo 

A simple integration by parts yields L e~ rt P(r > i) dt = (1 — E(e _rT ))/r for any nonnegative 

random variable r and r > 0. Therefore Problems flqlnfp and ( |qSup[ ) are equivalent to finding 

y{+) £ \> anc i y(-) g y respectively, such that 

(q±) sup ±E (e-rW) = ±E (W W±) )) , 

where ± denotes either + or — . Note also that if e g is an exponential random variable with 
E(e g ) = 1/g, independent of the filtration (J r t)t>o 5 then Problems ( |qlnf[ ) and ( |qSup[ ) are equivalent 
to minimising and maximising P(r(V) > e q ) over V G V, respectively. 
The following theorem holds. 



Theorem 1. A solution to Problem ( |q±[ ) is (for any q > 0) given by 

Remark 1. (i) Observe that by Theorem [IJ the mirror coupling V^ + ' = —B solves Prob- 
lem ( qlnfp and the synchronous coupling V^~' = +B is the solution to Problem flqSupP - 
(ii) Note that the solution depends neither on the parameters in @ nor on the discount rate q. 



3.1.1. Proof of Theorem{J\ Observe that, due to the symmetry in Problem ( |q±| ), we may assume 
without loss of generality that the starting points x, y in ©-(HJ satisfy (x,y) G D, where the set 
D C M 2 is given by 

(6) D :={(o,6); a > b > 0}. 

Fix q > and define the following function, closely related to the right-hand side in Problem ( |q±| ): 

(7) VP (±) (x,y) := E^ (e~^^) , (x,y) G D, 

where =b and =F simultaneously denote either + and — , or — and +. The proof of Theorem [1] is in 
two steps: we first establish sufficient conditions for a function \P : D — > M+ implying that ±\P is 
equal to the right-hand side in Problem ( |q=b[ ) (Lemmas [2] and [3]) , and then prove that fy^> in ([7]) 
satisfies these conditions (Lemma|3f. Throughout the paper we denote M + := [0, oo). 

For any measurable function ^ : D — > P + and Brownian motion V G V, consider the process 
U(y,V) = (U t (y,*))uilo,oo) defined by 

(8) U t (y,V) :=e-^(V))*(X tAT(v) ,Y tAT{v) (V)) 

(here and in the rest of the paper we denote s A t := min(s,t)). Then the following lemma (a 
suitable version of Bellman's principle) holds. 
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Lemma 2. Let ^f : D — > R + be a bounded continuous function satisfying fy(x,x) = 1 for all 
x > 0. //, for every (x, y) G D, the process ±U(V, \&) is a ¥ Xty -supermartingale for allV^V and 
U(^B,^) is a f x ,y -martingale, then V^ ' = =F-B solves Problem ( |q±| ), where ± and =f denote 
either + and —, or — and + . 

Proof. Since -X^y) = Y T ry\(V) P^-a.s. on the event {r(V) < oo} for any 7 G V, f is continuous 
and bounded, ^(x,x) = 1 holds for any x > and q > 0, the supermartingale property and the 
Dominated Convergence Theorem imply 

±E x>y (e-^W) = E Xiy (±U T{v) (V,V)I {T{v)<oo} ) < E X)V (±U (V,^)) = ±*(x,y), (x,y) G D, 

for all V G V (Ir.i denotes the indicator of the event {•}). Since t/(=F-B, \E f ) is a martingale, for 
V W = zs^B this inequality becomes an equality and the lemma follows. □ 

Our next task is to establish a verification lemma for Problem ( |q±| ). Let D° be the interior (in 
M?) of the set D defined in ([6]). For any twice differentiable function / G C 2 ' 2 (D°) we define the 
function /A '/ by the formula 

(9) [C (±) f) (x, y) := I axxf x + a 2 yf y + -cr 2 x 2 f xx + -cr 2 y 2 f yy =F o- 1 a 2 xyf xy - qf J (x, y), 

where (x,y) G D°, ± and =F denote either + and — , or — and +, and f x , f y , f xx , f yy and f xy 
denote the partial derivatives of /. For any function ^ : D — > R_j_, such that ^ G C 2,2 {D°) and 
Brownian motion VgV, the local martingale M(V, ^f) = (Mt(V, ^))te[o,oo)i given by 

(10) M t (V,*):= / e- qs {o- 1 X s * x (X s ,Y s (V))dB s + a 2 Y s (V)* y (X s ,Y s (V))dV s ), 

Jo 

is well-defined. 

Lemma 3. Assume the following hold: (I) $ : D —> M + is a bounded continuous function with 
V(x,x) = 1 for all x > 0; (II) * G C 2 ' 2 (£>°) and, in the interior D° , ^ xy < and C^H = 
(where ± denotes either + or — j; f/Z/J M(V, ^I 7 ) is a IF ' x ^-martingale for all (x, y) G D and V G V. 
TTien /or any (x, y) G D, V G V, £/ie process ±U(V,ty), defined in ([H]) ; is a F x ^ y - supermartingale 
and U(^B, \P) is a P X) „ -martingale. 

Proof. The definition of X and y(V) in ([3]) and Lemma LTOl in the Appendix imply d[X, Y(V)]t = 
C t aiX t a 2 Y t (V) dt, where C = (Ct)te[o,oo) is (J" t )-adapted and F(C t G [-1, 1]) = 1 for all t G [0, oo). 
Ito's lemma, the assumptions in Lemma [3] and definition (|8|) of U(V,^f) yield 



/•tAr(V) 
±l7 t (V,*) = ±*(ar,y)±M t (F,*) + / e- 9 Via 2 (l ± C s )X s y s (y)^ y (X s , Y s (y)) ds 

jo 



c*At(V) 
'0 

for all (x, y) G D and V G V. Since X, Y(V) and 1 ± C are non-negative processes and, by 
assumption @, we have o~io~2 > 0, the integrand in the representation of ±U(V, \&) is non-positive, 
making ±U(V, "if) a P^^-supermartingale. This representation, together with assumption (III), 
implies that U(^fB, \I / ) is a P^y-martingale. D 
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Recall that =p denotes either — or + and note the following equivalence: 

(11) P x<y (T(^B) = oo) = 1 for all (x,y) G D° ^^ q= = +, a 2 = a u a 2 < Oi- 

It is clear that under condition (jlip Theorem [1] holds. Lemmas [2] and [3] imply that in order 
to establish Theorem [1] in general, it is sufficient to prove that, when (|lip fails, the function 
Vl/W : D — )> K+ in ([7]) satisfies the assumptions of Lemma [3l More precisely, the following lemma 
holds. 

Lemma 4. Assumptions (I)-(III) of Lemma\3\ hold for the function VJ/W : D — > M + in (|7|), i//or 
some (x,y) £ D° we have Vx^t^B) = oo) < 1 (± and =F ore either + and — , or — and +j. 

Proof. Under the assumption of the lemma, the following representation holds: 

(12) M/(±)( x , y ) = g) fc± for (x,y) € A 
where 



A- 



-Ai/4 + ^W4) 2 + 2g/4, ifer±^0, 
a//i, if cr ± = 0, 

and a± and /i are defined in ©. Since, by assumption, the condition on the right-hand side 
in ([TT]) is not satisfied, the equality a± = implies /u > 0, making fc± a well-defined real number. 
Formula (|12p follows from the fact that r(=f B) has the same law as the first-passage time of the 
Brownian motion with drift (o~±Bt + yui)tg[o,oo) over the level log(x/y). The Laplace transform of 
this random time is given in [H p. 295] and amounts to the right-hand side of (|12p . 

Assumption (I) in Lemma [3] follows from (|12p . Furthermore it is clear that V^W £ C 2 ' 2 (D°). 
The formula in (|12|) and some simple calculations imply that for (x, y) £ D° the following holds: 

(13) *M(x,y) = - } ^*M(x,y), 9&(x,y) = ^V&(x,y), 

x y 

and 

k 2 / \ 

*w(s,y) = -* (±) (x,y) <0, (£ (±) tf (±) ) (x,y) = 0. 

y xy V / 

Hence assumption (II) of Lemma [3] is also satisfied. The equalities in (|13p and the definition 

in (|l(jp of the local martingale Af(V, \l/( ') imply that the integrands in the stochastic integrals 

are bounded processes and therefore square integrable. Hence M(V, VljW) is a P^^-martingale for 

all (x,y) £ D and V £V and assumption (III) of Lemma [3] also holds. □ 

3.2. Stationary problems. Consider the problems: find V mi £ V and V sup £ V such that 

(SInf) inf limsup ^ / P(r(F) > t) dt = limsup ^ / P(r(F inf ) > t) dt 

and 

(SSup) sup limsup ^ / P(r(y) > t) dt = limsup ^ f P(r(F sup ) > t) dt. 

A solution to these problems, independent of the values of the parameters of the geometric Brow- 
nian motions in ([3]), is given in the following proposition. Note in particular that, unlike in the 
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finite time horizon case, no new phenomena arise when the ergodic average criterion is used (i.e. 
the solution is completely analogous to the infinite time horizon case). 

Proposition 5. The Brownian motions V m ' = —B and V sup = B solve Problems (jSInfjl and QSSupP 
respectively. 

Proof. As in Section 13.1.11 we may assume that, due to symmetry, the starting points of X and 
Y(V) satisfy (x,y) £ D (see ©). By ([3]) and the definition of t(V) in Section [2] we have 

(14) t(V) = inf{t > 0; a 2 V t - a x B t + fit = log(x/y)}, 

where \jl is defined in ([5]) and the convention inf = oo is used. If x = y we have r(V) = for all 
V £ V and Proposition [5] follows. So we can assume (x,y) £ D° in the rest of the proof. 
We first analyse the case fi > 0. Since the Dominated Convergence Theorem implies 

limsup ^ / F(t(V) >t)dt= lim E((r(V)/T) A 1) = ¥(t(V) = oo), 

Problems (jSInfjl and ( |SSup[ ) are equivalent to finding V' ■ ' £ V such that 

(S±) inf ±P(t(V) = oo) = ±P(r(l/ (±) ) = oo). 

The strong law of large numbers for Brownian motion (e.g. [4, p. 53]), representation ()14p and 
log(x/y) > imply the equality F x ^ y (T(V) = oo) = for every V £ V and Proposition [5] follows. 

In the case /j<0we return to the formulation of Problems (jSInfj) and ( |SSup[ ) above. Observe 
that Theorem 0Jb]) yields the optimal couplings that minimise and maximise the probability 
P(t(V) > t) for every t > 0. Since the couplings are independent of t, they also minimise and 
maximise the stationary criteria in Problems (jSInfj) and QSSup[ ), which concludes the proof. □ 

Remark 2. The proof of Proposition \5\ relies in an obvious way on Theorem 0|b]) below. We 
would like to stress that there is no circularity in this argument since Proposition [5] is not used 
in Section [H Stationary problems are considered in Section [3] rather than later on in the paper, 
because the structure of the solution is the same as that of the infinite time horizon problems. 

4. Finite time horizon problems and the efficiency of the couplings 

4.1. Finite time horizon problems. Retain the setting and notation from Section [2j For any 
T > 0, consider the following problems: 

(T±) find V {±) £ V such that inf ±P(r(V) > T) = ±F(t(V^ ±) ) > T), 

where ± denotes either + or — . As in Section [31 we can reduce Problem (|T±[) to the case where 
diffusions in ([3]) start at (x,y) £ D, where D is given in ([6]). Define the set E := D x [0,T] and 
recall that a function F : E — > M + is the value function for Problem (|T±[) if: 

(15) F(x,y,t) = inf ±¥ x ^(t(V) > t) for all (x,y,t) £ E. 

Based on the results in Section [31 one might expect that ±<I > ( = ' = ', where 

(16) & ± \x,y,t):=F X)y (T(TB)>t), (x,y,t) £ E, 
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and db and =p simultaneously denote either + and — , or — and +, would be the value function for 
Problem (|Td=[) . In order to investigate this, we define the function A^- ' f for any / £ C 2,2,l {E°) 
(E° is the interior of E in R 3 ) by the formula 

[A^fj (x, y, t) := I aixf x + a 2 yf y + -o-\x 2 f xx + -cr 2 y 2 f yy =F o-\a 2 xyf xy - f t j (x, y, t), 

where (x,y,t) £ E°, ± and =F denote either + and — , or — and +, and f x ,f y ,ft, etc. denote 
the partial derivatives of /. For any sufficiently smooth function $ : E — > M_|_ and any Brownian 
motion V £ V, we define the local martingale N(V, <&) = (N t (V, &))te[o,T] by 

rtAr(V) 

(17) N t (V,$):= / (a 1 X s <f> x (X s ,Y s (y),T- S )dB s + a 2 Y s (V)$ y (X s ,Y s (V),T-s)dV s ). 

Jo 

The following proposition provides the key ingredient in the proof of Theorem [SJ 

Proposition 6. Let a bounded function <J> : E — > R + satisfy: (i) $(x, x,t) = for all x > and 
t £ [0,T], and $(x,y,0) = 1 for all (x,y) £ D° ; (ii) $ £ C 2 ' 2>1 (J5°) and, in the interior E° , i/ie 
equality A^'Q = holds, where ± denotes either + or — ; (wij i\T(F, $) is a W x ^ y -martingale for 
all (x,y) £ -D and V £ V. TTien t/ie following equivalence holds (± and =F denote either + and — , 
or — and +): 

&x y > on E° <s=>- V' ' = =pi? solves Problem (|Td=[) and d=<I> is its wa/ue function. 

Proof (=£•): The proof of this implication is analogous to that of Lemmas [2] (Bellman's principle) 
and [3] (submartingale property) in Section [3l The process ±U(V,$) = (±Ut(V, $))te[o,T\i 

(18) tft(V,$) := $(X tAriv) ,Y tAr(v) (y),T-t), 

is a P Xi y-submartingale for any V £ V and (x,y) £ .D (proof as in Lemma [3]). For any t £ [0,T], 
the boundary conditions in assumption (i) imply 

U t (V, $) = U T(K )(V; $) = P^-a.s. on {t > t(V)}. 

Hence, for any (x, y) £ D and V £ V, the submartingale property yields the inequality 

±P x , y (t(V) >T)= E Xiy (±U T (V, <S>)I {t{v)>t} ) = E x , y (±U T (V, $)) > ±E x>y U (V, $) = ±$(x, y, T). 

As in Lemma [21 this establishes the implication (note that, unlike Lemma [21 in this case we do 
not need, and in fact do not have, the continuity of <J> on E). 

(<=): Assume that there exists (xo,yo,7b) £ E°, such that $ X y(xo,yo,To) < 0, and that ±$ is the 
value function of Problem ( |Td=[ ), Bellman's principle implies that the process ±[/(V, <&), defined 
in (|18|) . is a ¥ x .^-submartingale for any V £ V and (x,y) £ D. Using our assumption, we now 
construct a Brownian motion V* ■ ' £ V, such that ±!7(V'' ',&) fails to be a P^^-submartingale 
(for any pair (x,y) £ D°), which will imply the proposition. 

The continuity of & xy implies that there exists r > 0, such that & xy is strictly negative on the 
set K 2 := H 2 x [T - 2r,T + 2r] C E° , where H 2 := [x - 2r,x + 2r] x [y - 2r,y + 2r]. Let 
.Hi := [xo — r, xo + r] x [yo — r, yo + r] and define the stopping times r x and r 2 by: 

r x (±) := inf{t £ [0,T]; (X,,F f (q=B)) £ H x }, r 2 (±) := inf{t £ [n,T}; (X t ,Y t (±B)) $ H 2 } 
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(where inf := T). Note that t[ ±] < r 2 (±) < T P Xjl/ -a.s. and Px^(t" 1 (±) < t 2 (±) ) > (there is a 
slight abuse of notation in the definition of r 2 as it is assumed that the process Y(±B), defined 
in fl«JJ), is driven by the Brownian motion ±B as indicated, but started at the random time r} 
and point Y (±)(q=B); ditto for X). 

r(±) 



Define the process yW = (V^ )tg[o,oo) ^y the following formula: 



v; 



(±) ._ 



„ , T W 1 (±) } ± V ) <^ ±) } T W^} 



dB s , 



where K/.i is the indicator of the event {•} and ± and =F denote either + and — , or — and +. 
Note that V"' • ' is an (Tt)-Brownian motion by Levy's characterisation theorem. Ito's formula on 
the stochastic interval [r{ , r 2 ] and assumptions (i)-(iii) in the proposition imply the following 
representation: 



E 



■'-,y 



±C/ (±) (y( ± ),$) + E 



x,y 



_(±) 



(±) 



2<j 1 u 2 X s Y s {V^)^ xy {X s ,Y s {V^),T - s)ds 



T 



(±) 



The event {r\ G (To — r, Tq + r),r(V r ( ± ^) > To + 2r} has strictly positive probability and the 
integrand under the conditional expectation is strictly negative on this event. We therefore find 

±C/ (± )(y (±) ,$)|T( ± )l <±[/ (± )0/ (±) ,$) onjr^ e (T -r,To + r),T(yW) >T + 2r} 

To ' T-i 7-1 



E* 



: v -a.s. This inequality contradicts the P x v -a.s. inequality 



z,2/ 



E.. 



±C7 (^(yW,^)!^ (±)1 > ±U (± )(F (±) ,$), 



which follows from the optional sampling theorem applied to the bounded P^y-submartingale 
U(y( ',§). This concludes the proof. D 

We will now apply Proposition [6] to study the question of whether ±$W, defined in (Tl6j) . is 
the value function for Problem (|T±[) . 

Lemma 7. Recall that [i and o± are given in ([5]) and assume a± ^ 0. Then, assumptions (i)-(iii) 
of Proposition^ hold for the function $W defined in (|16p . Furthermore, we have 



*gW*) = 21og(^)-y n f log(x/y)-/xf | + 



xy(|a ± |Vi) 3 



l°± 



Iv 7 * 



^yo-± Vy/ 



j.\W4 



A r 



log (x/y) - /it 



F± 



|V* 



for all (x,y) £ D° and t > 0, where N(-) is the standard normal distribution function and n(-) is 
its density. 

Proof. The explicit formula for the distribution of the running maximum of a Brownian motion 
with drift (see e.g. [3J p. 250]) yields the following representation of the function in (j 1 6 H : 

& ± \x,y,t)=h^ (log (x/y),t) for (x,y) e D, 



(19) 



where, for any z > and s > 0, we define 



h^\z,. 



N 



z — [IS 

W±\t/s 



oxi> I -^ 1 A" 



-Z — /US 

|cr±|v/i 
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Simple (but tedious) calculations using this representation yield the properties required in assump- 
tions (i)-(iii) of Proposition [6l For example, the martingale property of the process in (|17p follows 
from the fact that both x& x and y<&\ are bounded functions. The details are omitted. □ 

We are now ready to prove that the mirror (resp. synchronous) coupling of the driving Brownian 
motions in ([3]) is not necessarily optimal in Problem (T+) (resp. (T— )). In Theorem [HI we give a 
necessary and sufficient condition for the function ±$( ', defined (|16p . to be the value function 
for Problem flT±D , 

Theorem 8. Let ± and =F denote either + and —, or — and +, and recall that \i and o~± are 
given in ([5]) . Then the following holds for any time horizon and distinct starting points: 

(a) lf(i>0 and a± / 0, then V^ = =F-B does NOT solve Problem (1T±1) . 

(b) If /U < 0, then V^' = ^fB solves Problem ( |T±[) with the value function ±$ ( - ± - ) in (JTHJ) . 

Remark 3. (i) Note that under the assumptions of Theorem [FJJaj) , the mirror and synchronous 
couplings are suboptimal in Problems (T+) and (T— ) respectively. Furthermore, if ± = +, 
then o~± > and hence the optimality of the mirror coupling can fail even if the laws of X 
and Y(V) are equivalent (i.e. o\ = 02) for all V G V. 
(ii) In the case ^ > and a± = we have ± = — , o~\ = 0*2 and &(~\x,y, t) = I{t Ai <i og (x/y)} 
for all (x,y) G D°, t G [0,T] (recall (|14|) ). which implies that the synchronous coupling is 
suboptimal if and only if T > log(x/y)//i. 

Proof, (jaj) By Proposition [6] it suffices to show that for any fixed t > 0, there exists (x,y) G D° 
(see ([6]) for the definition of D) such that & xy (x, y, t) < 0. 

Define z := \og(x/y)/(\a±\\/t) > and a := fJ-Vt/\cr±\ > 0. Note that, since we are allowed 
to choose the point (x,y) G D° arbitrarily close to the diagonal half-line in the boundary of D, a 
Taylor expansion of order one of z *— > n(z—a) and z 1— > N{—z—a) around z = 0, the representation 
of &xy in Lemma [7] and the inequality 

(20) aN(-a) < n(-a) 

imply that <& X y (x,y,t) < for some (x,y) G D°. To check ([20]) . note that un{u) = —n\u) and 

/'OO poo 

aN(—a)= I an(u)du< / un(u) du = n(— a). 



(jb]) Assume first a± 7^ 0. Then the representation of Q xy in Lemma [7] and the assumption fi < 
imply & x ,y > on E° . Hence Proposition O yields the theorem. If a± = 0, we have ± = — , 
a 1 = a 2 and, by (HU), it follows ^~\x,y,t) = 1 for all (x,y) G D°,t G [0,T]. Hence -$H is the 
value function for Problem (T— ) and the theorem follows. □ 

4.2. Efficiency of the mirror and synchronous couplings. In this section we examine further 
the (lack of) optimality of the mirror and synchronous couplings characterised by the assumptions 
of Theorem [HJjaj). Since in this case the two couplings do not minimise and maximise (respectively) 
the coupling times of the geometric Brownian motions in §3§ over finite time horizons, but are 
nonetheless optimal both over the infinite time horizon (Section 13. ip and for to the stationary 
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criterion (Section I3.2p . it is natural to analyse whether the two couplings are efficient. A coupling 
V £ V is (exponentially) efficient (for some (x,y) £ D°) in Problem (1T±|) if the rate of the 
exponential decay of the tail of its coupling time is the same as the exponential decay of the value 
function F defined in (1151): 



lim - logP^ (t(V) >t)= lim - log (±F{x, y, t)) . 

t— s-oo t t— »oo t 

It is clear that if a coupling solves Problem (1T±|) for all time horizons T > and does not depend 
on T, then it is also efficient. Hence, by Theorem |8lfb|) , the mirror and synchronous couplings are 
efficient if \i < 0. However, the following statement holds. 

Proposition 9. If [X > and o~± 7^ 0, then the mirror and synchronous couplings are NOT 
efficient (for any (x,y) £ D°) in Problems (T+) and (T— ) respectively. 

Proof. The following bounds hold for the standard normal distribution function N(-), 

2 n(z) < N(z) < —z~ n(z) for any z < 0, where n = N' . 

i ~r z 

The first inequality follows from the identity /°°(1 + y~ 2 )e~ y > 2 dy = r _1 e _r > 2 for all r > 0, and 
the second is given in (|20p . These inequalities and the representation of 3H \x, y, t) in (|19|) imply 
that, for any (x,y) £ D°, there exist functions gi : (0, oo) — >■ (0, oo), i = 1,2, such that 

» W exp (- C° g '^> ~ "" 2 ) < *«(., ,, ( ) < „M exp (- "° g ^ t - ^ ) , 0, 

and lim^oo j log gi(t) = for i = 1,2 (the functions #j, i = 1, 2, also depend on a fixed starting 
point (x,y)). Since the logarithm is an increasing function, we find the limits 

1 1 u? 

lim -log(F x , y (r(TB) > t)) = lim - log($W(x,y, t)) = -f^-, 

i— ^oo t t— >-oo t 2,(7% 

where ± and =F denote either + and — , or — and +, which are independent of (x,y). 

Definition (0) and our assumption imply |ct+| > |<r_| > and hence /j,/(2aj t ,) < / u/(2a_). 
The mirror coupling is therefore not efficient for Problem (T+) since it has a strictly thicker 
exponential tail than the synchronous coupling. Likewise, the synchronous coupling is not efficient 
for Problem (T— ), which requires the thickest possible exponential tail among all couplings, since 
it has a thinner tail than the mirror coupling. This concludes the proof. □ 

Remark 4. It is the presence of the positive drift fi > that makes the mirror coupling suboptimal 
in Problem (T+) (see Theorem [8]). The proof of Proposition [9] suggests that if the drift fi is 
strictly positive and the time horizon T is large, it is in fact better (according to the exponential 
tail criterion) to use synchronous coupling. This naturally leads to the following conjecture: 

If jj, > and o± ^ 0, the synchronous (resp. mirror) coupling is efficient in Problem (T+) 
(resp. (T— )). 
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Appendix A. Family of Brownian motions on a filtered probability space 
Recall that V is defined in ([2]). See e.g. [6] for the proof of Lemma [101 

Lemma 10. For any Brownian motion V G V, there exists an {J- \) -Brownian motion W € V and 
a process C = (Ct)t>o, such that B and W are independent, C is progressively measurable with 
— 1 < Ct < 1 for all t > F-a.s., and the following representation holds: 

r-t ft 



V t = f C s dB s + [ T/T=C$dW, 
Jo Jo 



Remark 5. The proof of this lemma requires the existence of a Brownian motion B 1 - G V that 
is independent of B. If our probability space did not support such a Brownian motion, we could 
enlarge it, which would only increase the set V. Since the optimal Brownian motions in Theorems!]] 
andlHIfbj) are constructed from B alone, they would also have to be optimal in the original problem. 
Therefore we can assume that B exists. 
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